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Various theoretical uncertainties in the standard solar model and in the 
Mikheyev-Smirnov-Wolfenstein (MSW) analysis are discussed. It is shown 
that two methods of estimating the solar neutrino flux uncertainties are 
equivalent: (a) a simple parametrization of the uncertainties using the core 



method of Bahcall and Ulrich. In the MSW analysis, we emphasize proper 



^ i treatments of correlation of theoretical uncertainties between flux components 

H : 

5^ , and between different detectors, the Earth effect, and multiple solutions in a 

combined procedure. The MSW solutions for various standard and non- 
standard solar models are also shown. The MSW predictions of the global 
solutions for the future solar neutrino experiments are given, emphasizing the 
measurement of the energy spectrum and the day-night effect in Sudbury Neu- 
trino Observatory and Super-Kamiokande to distinguish the two solutions. 
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Typeset using REVT^ 



I. INTRODUCTION 



The solar neutrino experiments of Homestake (chlorine) [|T],3 , Kamiokande , and the 
gallium experiments of SAGE and GALLEX |^ show deficits of the neutrino flux from 
the Sun when compared to the standard solar model (SSM) predictions as summa- 
rized in Table |I[ Numerous theoretical proposals have been made to resolve the discrepancy 
between theory and experiment. Astrophysical solutions in general are, however, strongly 
disfavored by the data. As long as astrophysical processes cannot significantly distort the 
neutrino energy spectrum |T^ the lower observed Homestake rate relative to the Kamiokande 
rate excludes essentially all astrophysical explanations HTlHl^. Even with the Homestake 
experimental error tripled, the combined observations are in contradiction with the explic- 



itly constructed nonstandard solar models |11|. On the other hand, among many particle 



physics solutions, the Mikheyev-Smirnov-Wolfenstein (MSW) mechanism |jT5| gives an ex- 
cellent description of the data, and is taken as a strong hint of neutrino mass and mixings 
1|,|16|. (See also ^^^.) 

When constraining the parameter space from the experimental data in the MSW anal- 
ysis, it is necessary to include relevant theoretical uncertainties properly. In the SSM, the 
theoretical uncertainty in the initial flux quoted by Bahcall and Pinsonneault ^ is 14%, 
and is comparable to the experimental uncertainties of Homestake (10%) and Kamiokande 
(14%). In the Turck-Chieze-Lopes SSM the flux uncertainty is 25% and dominates 
the experimental uncertainties. The omission of the theoretical uncertainties underestimates 
the uncertainty of the MSW parameter space constrained from the experiments. 

Equally important, but often ignored, are the correlations among the theoretical uncer- 
tainties. Especially, a correct treatment of the ^B flux uncertainty is significant since it is the 
largest among the theoretical uncertainties and also strongly correlated from experiment to 
experiment. When one considers the combined fit, for example, it is not legitimate to allow a 
smaller ^B flux for Homestake and a larger ^B flux for Kamiokande. If the correlations were 
ignored in the MSW two-flavor oscillation analysis, one obtains a larger parameter space in 
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the large-angle solution and even finds a third allowed region around Am^ ~ 10~^ eV ^ and 
sin^ 26 ~ 0.7 at 90% CL. Moreover, the uncertainties are also correlated among different 
flux components. For instance, if the opacity were lower than the standard value (or equiv- 
alently the core temperature were lower), both the ''Be and '^B fluxes would be reduced. 
Since the MSW mechanism often affects each neutrino flux component differently according 
to the neutrino energy, the flux uncertainties and their correlations affect nontrivially the 
allowed parameter space of the combined observations. 

Two methods have been proposed to incorporate the flux uncertainties in the MSW anal- 
ysis. The flrst method utilizes the Monte Carlo SSMs constructed with randomly chosen 



input parameters distributed around the mean values [22]. Those 1000 SSMs were incor 



porated in the MSW analysis of Bahcall and Haxton |17|, and, recently, of Shi, Schramm 



and Bahcall 18 . Krauss, Gates, and White calculated the flux uncertainties from their own 



100,000 Monte Carlo SSMs ||2T|. [| This Monte Carlo method provides a rigorous estimation 



of the SSM uncertainties as long as their correlations are properly taken into account for 
different flux components and different experiments. 

The second method, which we have used in our previous analysis |T^|T^, parametrizes 



the SSM flux uncertainties with the central temperature and nuclear reaction cross sections, 
distinguishing the flux uncertainties due to purely astrophysical effects, such as the uncer- 
tainties from the heavy element abundance and other uncertainties in the opacity, from those 
due to the nuclear cross sections. The astrophysical uncertainties are parametrized by the 
uncertainty of the central temperature ATc, which is chosen to reproduce the uncertainties 
derived from the Monte Carlo estimation. This parametrization method has several advan- 
tages over the Monte Carlo method. The physical meaning of the uncertainties is clear, and 
it can be generalized to the nonstandard solar models that are in most cases parametrized 



^ The correlations among flux uncertainties are ignored in the analysis in Ref. [21|. Also the 



uncertainties used in Ref. |21| are larger than the estimation of Bahcall and Pinsonneault. 
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by the lower Tc |2^JTT|]. Also one can easily update the calculations by just changing the 
values of AT^ and nuclear cross section errors when a new solar model is introduced. It 
required only trivial changes to update our calculations from the Bahcall-Ulrich model to 
the Bahcall-Pinsonneault model, while the Monte Carlo SSMs are not yet available for the 
latter. 

This parametrization of the uncertainties, however, was questioned on the grounds that 
such a simplification can lead to errors in describing the nonlinear relations among the neu- 
trino fluxes, which are the output of solving the coupled partial differential equations of 
stellar structure with nontrivial matching conditions ||23| , p!^ . It was argued that a Monte 
Carlo study is necessary to estimate the uncertainties reliably. Later we will show numeri- 
cally that, in fact, the uncertainties obtained by he parametrization method are essentially 
identical to the Monte Carlo results. 

In this paper we discuss various technical but important issues concerning the theo- 
retical uncertainties in the solar neutrino data analysis. In Section we show that the 
parametrization method reproduces the SSM uncertainties obtained by the Monte Carlo 
method. The comparison is made for both the uncertainties of flux components and their 
correlations. We also compare the two methods for the uncertainties and their correlations 
of the predicted rates for different solar neutrino detectors. 

Once the equivalence of the two methods is established for the SSM, we compare them 



in the MSW analysis in Section |T|. The allowed regions calculated by both methods are 
shown. The effect of the theoretical uncertainties and their correlations are displayed for the 
global MSW analysis. Other theoretical issues are also considered in Section |T|. Analytic 
approximations for the MSW transitions by Petcov, Parke, and Pizzochero are compared and 
the associated uncertainties are discussed. Another issue involves the estimate of confidence 
level (C.L.) in the presence of multiple fit solutions. We emphasize that there are several 
possible definitions of the C.L. contours, leading to slightly different allowed regions, and 
give a statistical definition of the C.L. contours of the most conservative (and we believe the 
best) prescription. We also compare the allowed MSW parameter space for the SSM with 



increased uncertainties and the SSMs of different authors. 

We conclude Section |ITT| by presenting the current results of our MSW fits for transitions 
into both ordinary (z/^, Ur) and sterile neutrinos, incorporating the theoretical issues and 
uncertainties discussed above. These also include the regeneration in the Earth (Earth 
effect), which is important for some regions of the MSW parameters for both time-averaged 
data and day-night asymmetries p4| , p!6| . We find that the data is fit extremely well by the 
MSW effect for transitions into ordinary neutrinos for Am^ ~ 6 x 10~^eV^ and sin^ 2^ ~ 
7 X 10^^ (the nonadiabatic) solution, although there is a second (large-angle) solution with 
Am^ ~ 9 X 10^^ eV ^ and sin^ 29 ~ 0.6 which is marginally allowed at 90% CL. There is also 
a poorer but acceptable fit for transitions into sterile neutrinos in the nonadiabatic region. 

The MSW effect can be also considered in the nonstandard solar models, and the com- 



bined fits to explicitly constructed nonstandard solar models are shown in Section |V[ The 
results of the MSW fit using the core temperature and the *B flux each as a free param- 
eter are displayed. In Section |V], the prospects for the future solar neutrino experiments 
are considered. We discuss in detail the predictions for the energy spectrum measurement 
and the Earth effect in Sudbury Neutrino Observatory (SNO) and Super-Kamiokande to 
distinguish the two solutions obtained from the global analysis: a spectrum distortion is 
predicted for the nonadiabatic solution, while a characteristic day-night effect is expected 
for the large-angle solution. 

II. COMPARISON OF THE PARAMETRIZED UNCERTAINTIES AND THE 

MONTE CARLO RESULTS 

A. The Flux Uncertainties 



In our parametrization method [|T^,0, the basic assumption is that the SSM flux un- 
certainties are expressed by a few parameters that have physical meanings, and that those 
uncertainties are approximated by gaussian distributions. First we studied the distributions 
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of the 1000 Monte Carlo SSM fluxes calculated by Bahcall and Ulrich with randomly cho- 



sen input parameters |2^, and confirmed that they indeed have gaussian distributions. In 
Fig. |I], the histograms are the distribution of the Monte Carlo fluxes of the pp, ^Be, and ^B 
neutrinos; the solid lines are our fit with the gaussian form. The rate predictions for solar 
neutrino detectors were also studied and are displayed in Fig. 0, which are also fit with the 
gaussian form. In both cases, the fits are excellent. 

Secondly, we utilize the observation that the neutrino fluxes are well described by a power 



law in the central temperature 



Hpp) - Tc'-\ 0(Be) ~ and 0(B) ~ (1) 

where our units are such that Tq = 1 = 15.67 x 10^ K for the central value of the SSM. 
We express the astrophysical uncertainties of the major fluxes {pp, ^Be, and ^B) with the 
uncertainty of the central temperature times the exponent in the power law. We identify the 
main sources of the astrophysical uncertainties as the heavy element abundance and other 
uncertainties in the opacity. The other uncertainties, independent of such astrophysical ef- 
fects, are from the nuclear reaction cross sections, especially for the cross sections of p +^Be 
and ^He H-'^He. Those cross sections are expressed by S-factors 5*17 and 5*34, respectively. 
Our choice of the relevant cross sections is based on the fact that the ^B flux is directly 
proportional to Sn, and the ^Be and ^B fluxes are proportional to 6*34 The flux uncer- 
tainties and their correlations are expressed by those three parameters [] (ATc, ASn/Sn, 
and AS'34/S'34). Their contribution to each flux is displayed in Table ^ The magnitude of 



^ While the p +'^Be cross section has little effect on solar conditions other than the flux, other 
cross sections such as p + p (Su) and '^He +^He (/S'33) directly affect the energy generation and 
cannot be separated from the astrophysical uncertainties. In particular, the p + p reaction is the 
main source of the energy production, and a change of Sn leads to a change in Tc- We therefore 
consider that the effects of the uncertainties in and ^33 are included in ATq- 

^ We could extend this parametrization method by taking into account all nine SSM input pa- 
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the uncertainties of the major fluxes are the quadrature sum of those uncertainties: 

1/2 



A:=34,17 



i \2 



(2) 



~ T^'), and are the fractional 



where i = ppjBe,^B. Ui are the temperature exponents 
uncertainty in 0* from AS'i7/S'i7 and AS'34/S'34 listed in Table H for both the Bahcall-Ulrich 
and Bahcall-Pinsonneault model. AT^ is not explicitly given in the SSM calculation. We 
determine ATc so that the uncertainties defined above correctly reproduce the Bahcall- 
Ulrich or Bahcall-Pinsonneault flux uncertainties; the resulting ATc is consistent with the 
Tc distribution of the Monte Carlo SSMs displayed in Ref. We obtained ATc = 0.0057 
for the Bahcall-Ulrich model, and ATc = 0.0060 for the Bahcall-Pinsonneault model. (The 
Monte Carlo SSMs of the Bahcall-Pinsonneault model are not yet available.) 
The correlation matrix of the flux uncertainties is given by 



riirijy 



f A(f)' AcjP 



(3) 



fc=34,17 / \ 

where j = pp/Be,®B. 

The result of those uncertainties are compared to the Bahcall-Ulrich Monte Carlo SSMs, 
whose flux uncertainties are calculated by 



A0* 
0' 



N 



m=l 



1/2 



(4) 



rameters using the partial derivatives of the neutrino fluxes obtained from the Monte Carlo SSMs 
This should completely reproduce the Monte Carlo results. We will show, however, that our 
minimal choice of the parameters is sufficient to describe the SSM uncertainties. 
^ We similarly parametrized the minor fluxes [^^N, ^^O, "'^'^F, p + e + p {pep), and ^He + p {hep)] 
using Tc exponents n = 22, 28.5, 28.8, 2.8 and 4.5, respectively, and reproduce the amplitude of the 



uncertainties. However, the Tc exponents for these fluxes except hep are not obtained in Ref. |22|, 
and we do not reproduce the correlations of the pep and hep neutrinos with others properly. In the 
MSW calculations, the effect of the correlations among the minor fluxes is completely negligible. 



where (j)]^ are the ith fluxes of the Monte Carlo SSMs (m = 1,---,A^ = 1000); (f>^[= 
{J2m=i '^m) / ^] is mean value of the Monte Carlo fluxes. The correlation matrix is 
obtained by 

The results of the two calculations are compared in Table |T| for the magnitudes of 
the uncertainties and Table ^ for the correlations. The agreement of the magnitudes 
are excellent. The correlation matrices are also in good agreement especially for the pp-Be 
element. In Fig. ^, we display the distribution of Monte Carlo fluxes and our parametrization 
in 0(Be) — 0(B) plane; the agreement of the two methods are remarkable. 



B. The Rate Uncertainties 

The comparison of the two methods is also carried out for the predictions for different 
solar neutrino detectors. The fraction of contribution ff of the ith flux component to the 
dth detector {d =Kamiokande, CI, Ga) are listed in Table 0. The SSM uncertainty of the 
rate R'^ (in units of the central values of the SSM) for each detector is a quadratic sum of 
the flux uncertainties and the detector cross section uncertainties Aa'^/a'^, which is 0.033 
and 0.04 for the CI and Ga detector, respectively [p^ ; 

2 / . 2 



E (fl 

i=fiuxes \ 



Aa 



a" 



1/2 



(6) 



Those uncertainties are correlated by the ATc and si through the fluxes; the error matrix 
is 



ij=fluxes fc=34,17 



(7) 



where c,d = Kamiokande, CI, Ga. Here (and below) the correlation matrix Dd is related 
to the error matrix by 



D, 



^'^ AR^AR'^' 
8 



(8) 



Among the correlations, the flux uncertainty is most signiflcant because of its large am- 
plitude and strong correlation between the experiments, especially between the Kamiokande 
and CI rate. 

For the Monte Carlo method, the rate for the mth Monte Carlo SSM for the dth detector 
is given by 



(9) 



i=fluxes 

The rate uncertainties for the lith detector is 

I N 



m=l 



1/2 



(10) 



The error matrix element between the cth and dth detectors is 



N 



(11) 



m=l 



We compare the uncertainties and their correlations for the rates in Tables ^ (mag- 
nitudes) and VI I| (correlations). The parametrization method reproduces the Monte Carlo 
results with a remarkable accuracy. 



III. THE THEORETICAL UNCERTAINTIES IN THE MSW ANALYSIS 

A. The SSM Uncertainties 

We incorporate the SSM flux uncertainties described above in the MSW analysis using 
a method. The MSW calculations of the theoretical rate predictions for each experiment 
are described in Refs. [^,|T^]. (See also Refs. [P7|-P^.) The MSW rate for the dth detector 
{d = Kamiokande, CI, Ga) is 

<sw(sm'2^,A^')= E ftPf{sm'2e,Am'), (12) 

j=fluxes 

where Pf is the MSW survival probability of the fluxes {i represents nine flux components, pp, 
^Be(I), ^Be(II), *B, ^'^N, ^^O, ^''F, pep and hep) after integrating over the neutrino production 

9 



V,, = {ARt,,?+(^^K,^J + E (P^ff^)' (14) 



site and the neutrino energy including the detector cross sections; for Kamiokande the 
detector resolution and efficiency are also included when integrating over the recoil electron 
energy. The formula is also valid with the time-averaged Earth effect (but not with the 
Kamiokande II day-night data with six time bins). We calculate a value for each point 
in the sin^ 29 — Am^ parameter space; 

x'iAm^ sin^ 29) = ^ (i?,%, - i?^sw) (V-'U (Rt^t ' <sw), (13) 

c,d=Kam, CI, Ga 

where -Rexpt are the experimental values listed in Table |l[ is the 3x3 error matrix and 
its diagonal elements are the quadratic sum of the experimental uncertainties, the detector 
cross section uncertainties, and the SSM flux uncertainties: 

+ 

i=fluxes \ ^ / 

The off-diagonal elements describe the correlations of the flux uncertainties described by 
ATc and si; 

Va= E P^P,fU^[n'n^{ATcf+ E 44]- (15) 

i,j=fiuxes fc=34,17 

For the Monte Carlo SSMs, the is defined by a Monte Carlo average of the probability 
function: 

1 ^ 

exp(-xV2) = E exp(-X™/2) (16) 

m=l 

where Xm is the value calculated for the mth Monte Carlo SSM (m = 1, ■ ■ ■ , = 1000); 

( TDd _ Dd \2 
" d=Kam,Cl,Ga (^-^expt)^ + (-^m.MSW Ao'^ j O'^Y 

with the MSW predicted rate for the mth Monte Carlo: 

-^m.Msw ~ E ^^■^^~ir' (1^) 

i=fluxes T 
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We compare the Monte Carlo result and the parametrized method for both the Bahcall- 
Ulrich model and the Bahcall-Pinsonneault model in Fig. ^ when the Earth effect is ignored. 
The agreement between the Monte Carlo SSMs and the Bahcall-Ulrich model with our 
parametrization is excellent. 

It is important to include the possibility of v^. regeneration in the Earth [^,^. This 
affects the time-averaged rates and in addition, Kamiokande has searched for a day-night 
asymmetry by binning their data with respect to the angle between the nadir and the Sun. 
No asymmetry was observed. 



One can easily generalize Eqns ^ and 15 to include the Kamiokande day-night data 



point by expanding the error matrix to 9 x 9, representing the time-averaged rates 
of Homestake, gallium, Kamiokande III, and the six Kamiokande II day-night data points 
||16|| . We have scaled the normalized Kamiokande II day-night data taken from Ref. to 
the quoted Kamiokande II average value. Also we have added to each of the six bins the 
overall systematic uncertainty (15%) from the energy calibration, the angular resolution, 
and the event selection, which are factored out in the quoted normalized data. (We checked 
consistency by combining the six bins and reproducing the quoted average Kamiokande 
II rate.) The systematic uncertainty as well as the SSM flux uncertainties are properly 
correlated among the six bins and the other uncertainties. 

Fig. ^ shows the MSW allowed regions for the three SSMs when the Earth effect is in- 



^The 95% C.L. allowed regions in Figs. |, |, % % |, |T|, ||, |T|, ||, |T|, and |2| are defined by 
X^(sin^ 2^, Am'^) < Xmin + with Ax^ = 6.0 for both combined fits and individual experiment 
fits, which corresponds to Gaussian errors in two parameters. Improved definitions for Ax^ in the 
combined fits will be discussed later. For the individual experiments we use Ax^ = 6.0 (2 d.f.) 
instead of 3.9 (1 d.f.), since it is a more conservative estimate of the uncertainties and also easier 
to compare with the combined fits. (The latter would correspond to mapping a one-parameter 
confidence region onto a band in a two-parameter space.) 
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eluded. Again we conclude that the two methods yield essentially the same results, although 
the obtained large-angle region is slightly smaller in the parametrization method. 

We demonstrate the effect of the theoretical uncertainties by comparing the allowed MSW 
regions calculated without the flux and detector cross section uncertainties. Fig. |^(a) is 
the result without the theoretical uncertainties using the Bahcall-Pinsonneault SSM. When 
compared to Fig. ^(c), the nonadiabatic allowed region is noticeably smaller; there is no 
large-angle solution at 95% C.L. 

The omission of the correlation of uncertainties among the experiments can lead to a 
overestimation of the allowed parameter space in the large angle solution. To demonstrate 
the effect of the correlations, we have calculated the allowed regions without the uncertainty 
correlations, which is shown in Fig. ^(b). The correlations are significant in the large-angle 
region where the predicted Homestake rate is larger than its experimental central value, 
while the Kamiokande rate is smaller than its experimental central value, and the flux 
uncertainty does not enlarge the allowed region if the correlation between the two exper- 
iments are properly taken into account. Without the correlations, the allowed parameter 
space become larger. In Fig. ^(c) we display the allowed region ignoring both the correlations 
and the Earth effect, and also using the same experimental values used in Ref. ||2T]]. The 



large-angle solution stretches to Am^ as small as about 10 "^eV even at 90% C.L. Thus, 



the claim in Ref pT| that the allowed region is very large was in fact due to their neglect of 
the correlations. 

The allowed parameter space is significantly enlarged if the Turck-Chieze-Lopes SSM 
is used (Fig. ^, which predicts the smallest flux and assigns the largest uncertainties 
among the SSMs. We have also carried out the MSW calculations with doubled theoret- 
ical uncertainties [Fig. |^(a)], and with a tripled the Homestake experimental uncertainty 
[Fig. 1(b)]. 

The neutrino production distribution in the core and the electron density distributions 
in the Sun and Earth are potential sources of the uncertainties that are not included in 
the MSW calculations above. We have varied those quantities within reasonable range and 
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repeated the calculations, but no change is observed in the combined allowed parameter 
space. 



B. Analytic Approximations in the MSW Calculations 

Instead of solving the MSW differential equation numerically, we have calculated the 
neutrino survival probability with an analytic approximation proposed by Petcov that 
is obtained from an exact solution of the equation, assuming an exponential form of the 
electron density distribution in the Sun, which is a good approximation except for inside the 



core region (< 0.15 x Rq, where Rq is the solar radius). The Parke formula ||28|, another ap- 
proximation that is the simplest and assumes a linear electron density, yields essentially the 
same MSW plots as Petcov formula except for Am'^/E < 3 x 10~^eV^/MeV [see Fig. ||(b)], 
where the density variation is no longer well approximated with a linear function within the 
neutrino oscillation length. The adiabaticity proposed by Pizzochero assuming the exponen- 
tial electron density also gives a good approximation except for Am'^/E < 7x 10~^eV^/MeV 
[see Fig. §(c)], where the approximation fails because the oscillation length becomes larger 
than the density scale height at large angles. The three approximations are compared for 
the Kamiokande rate in Fig. differences are noticeable in the large angle with small Am^. 

One limitation of the Petcov formula (as well as the other two) is that it fails to de- 
scribe nonadiabatic level crossing when the neutrino production point is at or close to the 
resonant point. | As discussed by Krastev and Petcov [^, however, both of the global 
solutions shown in Fig. ^ are safe from this limitation. For the ''Be and neutrinos with 
energies relevant for the experiments, the resonance takes place at ~ 0.15Rq and ~ OARq, 
respectively, while most of those neutrinos are produced within O.IRq. The pp neutrinos 
are potentially dangerous since their energy corresponds to the resonance at or close to the 



^An analytic prescription when the neutrino production site overlaps with the resonance point is 
discussed by Haxton |31|. 
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center of the Sun, and there is a substantial overlap with the production site; however, in 
the allowed parameter space of the global fit, the MSW transitions of the pp are adiabatic 
and well described by the Petcov formula. We estimate that the uncertainty for the global 
fit due to the analytic approximation is < 1%. 



C. Confidence Level Definitions 



There are differences of the combined allowed MSW regions in the literature. Especially 



the large-angle regions shown in Refs. [ 18, 20, 21 1 are larger than our calculation, and even a 



third solution is allowed at 90% CL. for Am^ ~ 10~''eV^ and sin^ 2^ ~ 0.7. The possible 
sources of the difference are: treatment of the SSM flux uncertainties and their correlations, 
the treatment of the Earth effect, the MSW analytic approximations, and the experimental 
input data, some of which have been discussed above. The most significant difference, 
however, comes from the statistical definition of the confidence levels, rather than the details 
in MSW calculations. 

In the joint analysis above, we have used 

x\sm'2e,Am^)<xLn + ^x' (19) 

with Ax^ = 6.0 as the 95% confidence level (CL.) region in the 2-dimensional logsin^26' — 
logAm^ plane, where Xmm is the global minimum. This definition assumes a gaussian 
distribution of the probability density around the global minimum, which is only an ap- 
proximation in our case, particularly when we have multiple minima. Another definition 
used in the literature is to take the combined allowed regions simply as the overlap of the 
different experiments as shown in Fig. (See also Ref. pO|.) This definition can, however. 



overestimate the allowed regions. Consider, for example, the parameters that are marginally 
allowed at 95% CL. by two experiments; if we take the overlap as the combined fit, the pa- 
rameters are allowed. But the total x^ can be very large: x^ should be about 12 (= 2 x 6.0) 
at the edge of the two allowed regions and is allowed only at 99.8% CL. by the analysis. 



14 



Taking the overlap of the allowed regions of the different experiments displayed in Fig. ^ 
clearly overestimates the uncertainties when compared to the analysis shown in Fig. 
especially for the large-angle solution. Also the overlap procedure by definition ignores the 
correlations of the uncertainties between different experiments. 

Before improving the statistical definition of the uncertainties of the obtained parameters, 
we consider the goodness-of-fit for each MSW solution under the MSW hypothesis: if one 
of the MSW solutions is true, then how likely is it to obtain the observed (or a larger) x^l 
The goodness-of-fit is calculated from the minimum for each MSW solution. Without the 
Earth effect, the minimum of the nonadiabatic and large-angle solution are 0.5 and 4.9, 
respectively. For 1 degree of freedom (= 3 experiments - 2 parameters), the probabilities 
of obtaining the values equal to or larger than those values by chance are 48% and 
3%, respectively. That is, the hypothesis that the nonadiabatic solution is the true solution 
yields a good fit, while the large-angle solution hypothesis is possible statistically only at 3%. 
When the Earth effect and the Kamiokande II day-night result (6 data points) are included, 
there is a third x^ minimum in the large-angle region with sin^20 = 0.76 and Am^ = 
1.2 X 10"'' eV^, because of the regeneration of the electron neutrinos in the Earth at night 
in the pp and ^Be energy range. The x^ values are 3.1, 8.1, and 13.1 for the nonadiabatic, 
large-angle, and new large-angle solution, respectively. For 7 degrees of freedom (= 9 data 
- 2 parameters), the probabilities of getting x^ larger than those values are 88, 32, and 7%. 
The fits for the first two solutions are reasonable, while the third fit is somewhat poor but 
not excluded. The small x^ contribution from the Kamiokande day-night data is responsible 
for the improvement of the large-angle, large Am? solution. The results are summarized in 
Table [VII1| and |X| and the corresponding results for sterile neutrinos in Tables ^ and 



The nonadiabatic solution gives a better fit than the large-angle solutions, either with 
or without the Earth effect, and this tendency can be quantified in another way: suppose 
the two-flavor MSW is true and the probability density of finding the true parameter is 
distributed throughout the log sin^ 26 — log Am^ plane, including all the x^ minima. That 
is, we assume that the probability distribution is 
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P{sm'^2e,Am'^) = exp[-x^(sin^ 2^, Am2)/2], (20) 

where is chosen so that the total probabihty is unity. Then what is the relative probabili- 
ties of finding the true parameters in the diff'erent regions? We approximate the probability 
distributions as an overlap of Gaussian distributions, each corresponding to a different MSW 
region. Then the relative probability of finding the true parameters in the ith region [i = 
nonadiabatic, large-angle solutions) is calculated as 

pi^^ 27ra>;,v/r^exp(-xf,^i./2) ^^^^ 

re ative 2TTaiaLV'^ - P^'^ exp(-Xj^,mm/2) ' 

where as and am are the standard deviations of logsin^26' and logAm^, respectively, and 
is the correlation parameter; j runs over each allowed region with the minimum Xjmin- 
Without the Earth effect, the relative probabilities of finding the true parameters in the 
nonadiabatic and large-angle region are 89% and 11%, respectively. When the Earth effect 
and the Kamiokande II day-night data are included, the probabilities for the nonadiabatic 
and two large-angle solutions are 94.9%, 4.6%, and 0.5%, respectively. The nonadiabatic 
solution is strongly favored. The results are summarized in Table |V111| and ^X. 



Finally, using the probability density approximated with multiple Gaussian distribu- 
tions, we improve the confidence level definition of the MSW parameter uncertainties. The 
confidence level P and Ax^ are related by 



Ej 2n aiain^X - (P^ exp(-x^,nim/2) 



where 

^'^amm'^^^ r 1 



V.{Ax') = P:"^ " exp[-i(x^,_ + x'- xLin)] dx' (23) 

with Xomin = '^i'^l X^min N = all local minima}. Without the Earth effect, the 90, 95, 
and 99% C.L. correspond to Ax^ = 5.9, 7.3, and 10.6, respectively. With the Earth 
effect, Ax^ = 5.5, 7.0, and 10.9. (These should be correspond with the values Ax^ = 
4.6, 6.0, and 9.2 for a single Gaussian distribution.) We believe that the new prescription 
is more reliable because it takes into account the existence of multiple minima and our lack 
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of a priori knowledge of which is the true solution. It is also a more conservative estimate 
of the uncertainties. The allowed regions with the improved CL. definition are shown in 

Fig. in. 



D. Results of the Global MSW Analysis 



The results of the global MSW analysis for two-flavor oscillations are displayed in Fig. ^ 



and the best fit parameters are listed in Tables These include the SSM uncertainties and 



their correlations, the Earth effect and the Kamiokande II day-night data, the improved sta- 
tistical definition of confidence levels, and the updated experimental input data. [| Between 
the two allowed regions at 90% C.L., the nonadiabatic solution with 

Am^ ~ 6 X lO^^eV^ and sin^ 2^ ~ 7 x 10~^ (24) 

gives an excellent description of the observations, while the second (large-angle) solution 
with 

Am^ ~ 9 X 10"'^ eV ^ and sin^ 29 ~ 0.6 (25) 

is marginally allowed at 90% CL. In the second region the Earth effect and the uncertainty 
correlations between the experiments are significant: the regeneration of z/g in the Earth 
during the night distorts and enlarges the allowed region while the absence of the day-night 



asymmetry in the Kamiokande II data excludes a wide parameter space |T6| , |24| , p3| . The 
omission of the correlation would result in overestimating the uncertainties. There is also 
a third solution at large-angle and low Am^ that is allowed at 99% CL. when the Earth 
effect, which in this region is significant for the pp and ^Be fluxes, is included. 



^There are additional Kamiokande II data of the energy spectrum of recoil electrons Q. The spec- 
trum slightly disfavors a part of the adiabatic region (Am^ ~ 10""^ eV^ and sin 2^ ~ 6 X 10 ), 
but the effect in the combined fit is insignificant. The spectrum information is not statistically 
independent and is therefore not included in our analysis. 
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One can consider the MSW oscillation to sterile neutrinos instead of z/^ or Vr- In that 
case the survival probability of electron neutrinos depends both on the electron and neutron 
density in the Sun |1^J2^J3^. A more important difference is the lack of a neutral current 
contribution in the Kamiokande experiment, which would amount to about 15% of the total 
signal for flavor oscillations, requiring a larger v^, survival probability in the sterile neutrino 
case. As a result the MSW effect now cannot completely resolve the larger Kamiokande 
rate relative to that of Homestake, and the fit becomes poorer, especially for the large-angle 
(with large Am^) solution. There, Xmin is 21.0 for 7 degrees of freedom including the Earth 
effect, and no allowed region exists even at 99% C.L. (see Fig. |13|, 0, Tables 0, and |X| .) 



The large-angle region for sterile neutrinos is independently excluded by the bound on the 
number of neutrino species in big-bang nucleosynthesis 



IV. THE MSW EFFECT IN NONSTANDARD SOLAR MODELS 

We have so far considered the MSW effect within the SSM uncertainties. The SSM is, 
however, still a theory that needs to be calibrated. Although the nonstandard solar models 



alone cannot explain the solar neutrino data [Tl||, it is possible that the MSW effect takes 



place while the SSM is incorrect. We have examined constraints on the MSW parameter 
space with four different nonstandard solar models that are explicitly calculated and have 
neutrino flux predictions. Each of these models is ad hoc and assumes nonstandard input 
parameters that are grossly different from those in the SSMs, addressing possibilities of 
our ignorance of the astrophysical quantities such as the opacity or the primordial element 
abundances, or of the nuclear reaction cross sections that have never been measured at 
energies equivalent to the solar temperature. Two of them were constructed in attempts to 
explain the solar neutrino deficit, and predict smaller '^B and ''Be fluxes than the SSMs: the 
low opacity model by Dearborn ||36|, in which the opacity is reduced by 20% in the region 
where the solar temperature is above 5 x 10^ K, and the high Su model with the p + p cross 
section increased by 30% The allowed regions are generally enlarged compared to the 
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SSM case, and even the adiabatic branch of (Am^ ~ 10~^ eV^ and sin^ 26 ~ 10~^ — 1) 
is allowed. The MSW parameters for the two models are shown in Fig. |15|(a) and (b). Q 
The two other nonstandard solar models predict larger neutrino fluxes: the high Y model 



2^1 that assumes a larger primordial helium abundance in the solar interior, considered to 
solve problems with helioseismology data, and the maximum rate model in which 6*33 is 
artificially set to zero, yielding the largest prediction for the gallium rate (303 SNU). With 
the larger fluxes the combined allowed regions generally move inward in the MSW triangle, 
as displayed in Fig. |l6|(a) and (b). 

The solar models with a nonstandard opacity value or nonstandard heavy element abun- 



dance are generally parametrized by a nonstandard central temperatures (Tc) [p!T| -[T3[|, and 
the neutrino fluxes are described by the power laws of Eqn. (^. One can constrain the 
MSW parameter space from the combined data for these generic nonstandard solar models 



by allowing Tc to change freely. The allowed MSW regions are displayed in Fig. |T^. The 
existing observations constrain 

Tc = 1.02 ±0.02 (la), (26) 

with x^jj^ = 2.1 for 6 degrees of freedom (d.f.); the Tc obtained by observations allowing the 
MSW effect is consistent with the SSM (Tc = 1 ± 0.006). We note that without the MSW 



effect there is no consistent Tc to describe the observations |]13| , p!2| , p | . 

Another major source of uncertainty in the SSM is the p+^Be cross section (S'17), which 
is poorly measured and has the largest uncertainty among the SSM input parameters. S'17 
is directly proportional to the ^B flux and is independent of the astrophysical uncertainties 
parametrized by Tc- We have carried out MSW fits assuming two extreme values of Sn] 
a 30% reduction and a 50% increase. The results are shown in Fig. [l^(a) and (b). It is 
also reasonable to consider the ®B flux as a free parameter, given the large uncertainties 



°No flux uncertainties are given in those models. We used the standard values of the uncertainties 
in Ref. [g. 
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from both astrophysics and nuclear physics. When the flux is used as a completely free 
parameter in the MSW flt, the experiments constrain 

</.(B)/0(B)ssM = 1.43 + 0.65 - 0.42 (la) (27) 

with Xmin = 2.1 for 6 d.f., and the allowed regions are displayed in Fig. |TB|(c). The larger 
values of 0(B) and Tc are preferred since a larger '^B flux can reduce the relative difference 
in the survival probabilities for Homestake and Kamiokande, giving more freedom in the 
MSW parameter constraint. 



V. MSW PREDICTIONS FOR FUTURE EXPERIMENTS 

One can check the consistency between the current observations and the MSW theory by 
predicting the gallium rate from the combined chlorine and Kamiokande data pSf . The pa- 



rameter space allowed by the two experiments is shown in Fig. |T^(a), along with the survival 
probability contours for the gallium experiment. Assuming the Bahcall-Pinsonneault SSM, 
the combined Homestake and Kamiokande observations predict the gallium rate to be < 100 
SNU; the MSW prediction is consistent with the current result of SAGE and GALLEX al- 
beit the large theoretical uncertainty. By reducing the statistical uncertainties, we estimate 



the allowed parameter space expected in the near future. Shown in Fig. |T9|(b) is the region 
allowed by the combined observations when the gallium statistical uncertainty is reduced 
by a factor 1 / y^, which is equivalent to the data set of the gallium experiment by the end 
of 1994; the central value is assumed to be the same as current data. Although the nona- 
diabatic region shrinks somewhat and the large-angle solution is no longer allowed at 95% 
C.L., no significant change is expected unless the experimental values change drastically. 

Although the current observations are consistent with the two-flavor MSW hypothesis 
and strongly disfavor the astrophysical solutions, a smoking gun evidence for a nonstan- 
dard neutrino physics is still awaited. Theoretical questions yet to be answered are: (a) 
distinguishing astrophysical solutions and particle physics solutions (e.g., the MSW effect), 
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(b) calibrating solar models, and (c) distinguishing the MSW nonadiabatic and large-angle 
solution and extracting Am^ and sin^ 26 from the data. 

For the next-generation solar neutrino experiments, such as Sudbury Neutrino Obser- 
vatory (SNO) H, Super-Kamiokande |0|, BOREXINO and ICARUS [0, the MSW 



mechanism can yield robust predictions and is a verifiable hypothesis. The measurement of 
the charged to neutral current ratio (CC/NC) in SNO and the measurements of the neutrino 
energy spectrum and the day-night effect in SNO and Super-Kamiokande should be able to 
answer the questions addressed above. 

The measurement of the neutral current in SNO is insensitive to flavor oscillations and 
is a direct measurement of the *B flux. One can calibrate the core temperature of the Sun 
at the 1% level given the power law dependence of the flux (0(B) ~ T^) and allowing a 
20% uncertainty from the experimental error and the Si-j uncertainty. 

The measurement of CC/NC in SNO would be the most direct test of the MSW hypoth- 
esis. The up-to-date global analysis predicts the ratio to be 

{0.2 — 0.6 (nonadiabatic solution) 
(28) 
0.2 — 0.3 (large-angle solution) 

compared to the ratio expected if no oscillations occur. This measurement is, however, 
insensitive to oscillations between electron neutrinos and sterile neutrinos. 

The measurements of the charged current spectrum in SNO and the recoil electron spec- 
trum in Super-Kamiokande are important for distinguishing the two solutions suggested by 
the global analysis of current observations. Should the nonadiabatic solution, which gives 
the better flt for the data, be the case, we expect a depletion of electron neutrinos in the 
lower end of the observed spectrum, while little distortion is expected in the large-angle 
solution (or for astrophysical solutions). In Fig. we compare the spectra of the two 
solutions, one in the best flt solution in the nonadiabatic region, and one in the large-angle 
region. The large-angle spectrum is normalized to the nonadiabatic spectrum above the 
threshold. The error bars indicate the statistical uncertainties assuming 6,000 events, equiv- 
alent to two years of operation. We have included in the calculation the charged current 
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cross section [^3| , p| and the detector energy resolution In Fig. |2T|, the spectra of the 
two solutions are shown for 16,000 events (two year operation) in Super-Kamiokande. The 
spectrum shape expected for oscillations to sterile neutrinos are almost identical to these. 

Another MSW prediction that helps to distinguish the two solutions is the day-night 
effect due to matter oscillations in the Earth during the night. The effect shows up not only 
in day-night differences of signals, but also in the time dependence during the night and in 
seasonal variations due to the obliquity of the Earth. The large-angle solutions are quite 
sensitive to the Earth effect, but the nonadiabatic solution is insensitive at the observable 
level. f\ The expected rates at night along with the day-time rate are shown in Fig. ^ and 



23| for SNO and Super-Kamiokande, respectively. The night rate is divided into six bins 
according to the angle between the Sun and the nadir at the detectors, corresponding to 
bin neutrinos with different path lengths in the Earth. The effect yields such a noticeable 
variation in the large-angle region that an updated Kamiokande day-night data might be 
enough to confirm or rule out the large-angle solution. The prediction for Kamiokande with 
error bars equivalent to 200 events is shown in Fig. E^. 



The spectrum distortions and the day-night effect are particularly important for oscil- 
lations to sterile neutrinos, for which the charged to neutral current ratio is unchanged by 
the oscillations and is therefore the same as for astrophysical solutions. The distortion of 
the spectrum signifies the nonadiabatic solution, while the day-night effect indicates the 
large-angle solution, which is already excluded at 99% CL. by the existing data. 

When the SSM and two-flavor oscillations are assumed, one can predict the rate in 



those experiments. In Fig. 25, the combined allowed parameter space is displayed with the 



possibility of detecting the day-night effect for the nonadiabatic solution was discussed in 
Ref. |^5|. While it is true that the instantaneous enhancement of the night-time signal can be 
~ 20% with respect to the day-time signal, it is hardly measm'able when the signals are averaged 
over time bins and the statistical uncertainties are taken into account. See Fig. ^(a) and ^(a). 
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survival probability contours of SNO (the charged current mode), Super-Kamiokande, and 
BOREXINO. Those high-counting experiments should be able to constrain the parameter 
space precisely and check the consistency of the MSW predictions. 



VI. CONCLUSIONS 

Various theoretical uncertainties in the global MSW analysis have been discussed in 
detail. It was shown that our parametrized SSM uncertainties yield essentially the same 
result as the Monte Carlo estimation by Bahcall and Ulrich. The direct comparison was 
made in the flux uncertainties and their correlations, in the rate uncertainties and their 
correlations, and in the MSW calculations. The different MSW approximations of Petcov, 
Parke, and Pizzochero were compared. Various confidence level definitions are discussed and 
a careful statistical treatment in the global fit was emphasized. 

There are two MSW global solutions at 90% C.L., one in the nonadiabatic region and 
the other in the large-angle, the former solution giving a considerably better fit. The proper 
treatment of the Earth effect and the Kamiokande day-night data is significant in the large- 
angle region. When the Earth effect for the ^Be and pp neutrinos are included, there is 
a third solution at 99% CL in the large-angle, small Am^ region. For the oscillations 



to sterile neutrinos, the solution is limited to the nonadiabatic at 90% CL. In Fig. ^ the 
MSW solutions are displayed with other observational hints of neutrino mass: the oscillation 
interpretation of the atmospheric neutrinos |^ and the cold plus hot dark matter scenario 
to interpret the Cosmic Background Explorer (COBE) microwave background anisotropy 
measurement and the large-scale structure observations E^. The see-saw predictions for 



neutrino mass and mixing are also displayed in Fig. 26: t'^ oscillations in the SOio 



grand unified theory (GUT) with an intermediate-scale breaking |]12|,B8|, z/g i^t oscillations 



in the supersymmetric SOio GUT and Ve ^ oscillations in the superstring- 



inspired model with nonrenormalizable operators |^9|. In the SOiq GUT, the i>r mass is 



expected to be in a range relevant to the cosmological hot dark matter. Assuming the 
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Bahcall-Pinsonneault SSM, the MSW solutions for Am^ are generally in agreement with 
the theoretical expectations, especially with the v^. *^ oscillations in the supersymmetric 
SOio model. However, the mixing angles are not consistent with the expectation of the 
simplest versions of the models that the lepton mixing angles are similar to the corresponding 
quark mixing angles [O]. (For the string-inspired model of Ref. p3 there is no compelling 
prediction for the mixing angles.) 

The global MSW was also carried out for various standard and nonstandard solar models; 
with the Turck-Chieze-Lopes SSM or nonstandard solar models a wider range of the param- 
eter space is possible. When the core temperature and the flux each was used as a free fit- 
ting parameter, the data constrained Tq = 1.02±0.02 and 0(B)/0(B)ssm = 1.43 + 0.65 — 0.42 
at la, respectively. 

The predictions of the MSW solutions assuming the Bahcall-Pinsonneault SSM were 
discussed in detail for SNO and Super-Kamiokande. We expect for flavor oscillations the 
charged to neutral current ratio in SNO to be 0.2 - 0.6 and 0.2 - 0.3 of the SSM prediction 
for the nonadiabatic and large-angle solution, respectively. The nonadiabatic solution yields 
spectrum distortions measurable in SNO and Super-Kamiokande, while the large-angle solu- 
tion predicts characteristic day-night differences. The spectrum distortion and the day-night 
effect are independent of solar models, and are particularly important for the oscillations 
to sterile neutrinos because the absence of the neutral current in SNO prevents one from 
distinguishing neutrino oscillations from astrophysical solutions in the charged to neutral 
current ratio measurement. 
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TABLES 

TABLE L The standard solar model predictions of Bahcall and Pinsonneault (BP SSM) and 
of Turck-Chieze and Lopes |^ (TCL SSM), along with the results of the solar neutrino experiments. 
The gallium experiment is the combined result of SAGE and GALLEX I and IL 



Kamiokande ^ 

Homestake ^ (SNU) 

SAGE ^ & GALLEX ^ (SNU) 



BP SSM TCL SSM 

1 ± 0.14 0.77 ± 0.19 

8±1 6.4 ± 1.4 

131.5 122.5 ± 7 



Experiments 
0.51 ± 0.07 BP SSM 
2.32 ± 0.23 (0.29 ± 0.03 BP SSM) 
81 ± 13 (0.62 ± 0.10 BP SSM) 



^The result of the combined data of 1040 days of Kamiokande II [0.47 ± 0.05 (stat) ± 0.06 (sys) 
BP SSM] and 514.5 days of Kamiokande III [0.57 ± 0.06 (stat) ± 0.06 (sys) BP SSM] @. 
'"The result of Run 10 to 124 (through May, 1993) §. 

'^The prehminary result of SAGE I (from January, 1990 through May, 1992) is 70 ± 19 (stat) ± 10 
(sys) SNU [§. 

'^The combined result of GALLEX I and II (including 21 runs through April, 1993) is 87 ± 14 
(stat) ± 7 (sys) SNU 0. 
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TABLE II. The parameters of the SSM flux uncertainties, n is the exponent in the Tc power 
law, and Sk = AS^/Sk (k = 34, 17) are the fractional flux uncertainties due to the nuclear reaction 



cross section '^He+^He and p+^Be. We quote from Ref. |22,||. We determine ATq to be 0.0057 
for Bahcall-Ulrich SSM and 0.0060 for the Bahcall-Pinsonneault SSM. 



n S34 Sir 

Bahcall-Uhich SSM {ATc = 0.0057) 

pp -1.2 

^Be 8 0.02 

^B 18 0.02 0.07 
Bahcall-Pinsonneault SSM {ATc = 0.0060) 

pp -1.2 

^Be 8 0.032 

^B 18 0.032 0.093 
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TABLE III. The magnitudes of flux uncertainties {A(j)/(p at la) quoted from the Bahcah-Uhich 



SSM |22|, of the Bahcall-Ulrich Monte Carlo SSMs (Gaussian fit), and of the parametrized Bah- 
call-Ulrich SSM using the central temperature and the nuclear reaction cross sections. Also listed 
are the uncertainties of the Bahcall-Pinsonneault SSM and of its parametrized fluxes. (The Monte 
Carlo study of the Bahcall-Pinsonneault model is not available.) 



PP 



7 



Be 



^B 



Bahcall-Ulrich SSM 
Bahcall-Ulrich SSM (Monte Carlo) 
Parametrized (ATc = 0.0057) 



0.0059 
0.0067 
0.0069 



0.050 

0.05 

0.05 



0.12 
0.17 
0.13 



Bahcall-Pinsonneault SSM 
Parametrized (ATc = 0.0060) 



0.0067 
0.007 



0.06 
0.06 



0.14 
0.15 



TABLE IV. The correlation matrices of flux uncertainties obtained from the Bahcall-Ulrich 
Monte Carlo SSMs and the parameterization method. The agreement between the two methods is 
good, especially for the pp ~ ^Be element. 



PP 



'Be 



Bahcall-Uhich SSM (Monte Carlo) 

PP 
8B 
^Be 



1 

-0.73 
-0.92 



1 

0.74 



Parametrized with ATc and As 



PP 
8B 
■^Be 



1 

-0.81 
-0.92 



1 

0.80 
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TABLE V. The fractional contribution from flux components for the predicted SSM rates for 
the different solar neutrino detectors. The total SSM value is normalized to one. 



Bahcall-Ulrich SSM 

pp 0.536 

^Be (I) 0.009 

^Be (II) 0.139 0.251 

^B 1 0.772 0.106 

pep 0.025 0.023 

0.013 0.029 

^^O 0.038 0.046 

Bahcall-Pinsonneault SSM 

pp 0.538 

^Be (I) 0.009 

^Be (II) 0.150 0.264 

^B 1 0.775 0.105 

pep 0.025 0.024 

0.013 0.023 

0.038 0.037 
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TABLE VI. The comparison of the magnitudes of rate uncertainties for the Bahcall-Uhich 
SSM, Monte Carlo SSMs, and the parametrized SSM. Also listed are the uncertainties in the 
Bahcall-Pinsonneault SSM and its parametrized SSM. The detector cross section uncertainties are 
included in the chlorine (CI) and gallium (Ga) uncertainties. 





Kamiokande 


CI 




Ga 


Bahcall-Ulrich SSM 


0.12 


0.11 




+0.05 - 0.04 


Monte Carlo 


0.12 


0.11 




0.05 


Parametrized (ATc = 0.0057) 


0.13 


0.11 




0.05 


Bahcall-Pinsonneault SSM 


0.14 


0.13 




0.05 


Parametrized {ATc = 0.0060) 


0.15 


0.13 




0.05 


TABLE VII. The comparison 


of the correlation matrices of rate uncertainties for the Bah- 


call-Ulrich Monte Carlo SSMs and the parametrized SSM. The agreement 


is excellent. 




Kamiokande 




CI 


Ga 


Bahcall-Ulrich SSM (Monte Carlo) 










Kamiokande 


1 








CI 


-0.997 




1 




Ga 


-0.92 




0.95 


1 


Parametrized with ATc and As 










Kamiokande 


1 








CI 


-0.996 




1 




Ga 


-0.90 




0.94 


1 
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TABLE VIII. The best fit values of the MSW parameters when the Earth effect is ignored. 
P is the goodness-of-fit, i.e., the probabihty of obtaining by chance a equal to or larger than 
the obtained x^- -^relative is the relative probability between the different solutions when the 
probability distribution is gaussian for each solution. The improved definition of 90, 95, and 99% 



C.L. correspond to Ax^ = 5.9, 7.3, and 10.6 [see Fig. |Tl|(a)]. 







Nonadiabatic 


Large Angle I 


sin^ 29 
Am2 (eV2) 




6.4 X 10-3 
6.3 X 10"^ 


0.71 

1.6 X 10"^ 


(1 d.f.) 

P(%) 

-^relative (%) 




0.5 
48 
89 


4.9 
3 
11 


TABLE IX. Same as Table VIII except that the Earth effect and the Kamiokande II day-night 
effect are included. The improved definition of 90, 95, and 99% C.L. correspond to Ax^ = 5.5, 7.0, 
and 10.9 [see Fig. |ll|(b)]. 




Nonadiabatic 


Large Angle I 


Large Angle II 


sin^ 29 
Am2 (eV2) 


6.5 X 10^3 
6.1 X 10~6 


0.62 
9.4 X 10-6 


0.76 
1.2 X 10"^ 


(7 d.f.) 

P{%) 

-^relative (%) 


3.1 
88 
94.9 


8.1 
32 
4.6 


13.1 
7 

0.5 
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TABLE X. Same as Table |VIII except that the oscillations are for sterile neutrinos. The Earth 
effect and the Kamiokande II day-night effect are not included. The improved definition of 90, 95, 
and 99% C.L. correspond to Ax^ = 4.8, 6.5, and 10.6. 



Nonadiabatic Large Angle I 



sin^ 29 




7.4 X 10^3 




0.85 


Am2 (eV2) 




4.7 X 10"^ 




9.6 X 10^6 


(1 d.f.) 




2.8 




11.7 


P(%) 




9.4 




0.06 


-^relative (%) 




98.9 




1.1 


TABLE XL 


Same as Table ^ except that the Earth effect and the Kamiokande II day-night 


effect are included. The improved definition of 90, 95, 


and 99% C.L. correspond to Ax^ = 4.6, 6.0, 


and 9.1. 












Nonadiabatic 




Large Angle I 


Large Angle II 


sin^ 29 


7.0 X 10^3 




0.77 


0.60 


^rr? (eV2) 


4.5 X 10~6 




6.7 X 10-6 


6.9 X 10"^ 


X' (7 d.f.) 


7.0 




21.0 


24.2 


P(%) 


43 




0.4 


0.1 


^relative (%) 


99.59 




0.03 


0.01 
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FIGURES 

FIG. 1. The distributions of the (a) pp, (b) ^Be, and (c) flux of the 1000 Monte Carlo SSMs 
by Bahcall-Ulrich (histograms) are fit with the parametrized method (sohd curves) that assumes 
Gaussian distributions of the central temperatures and the nuclear reaction cross sections around 
their central values. 

FIG. 2. The distributions of the experimental rates obtained from Bahcall-Ulrich Monte Carlo 
fluxes for the (a) Kamiokande, (b) chlorine, and (c) gallium experiments (histograms). They are 
compared with the rate distributions obtained from the parametrized fluxes (solid curve). In both 
cases the detector cross section uncertainties are not included. 

FIG. 3. The distributions of the ^Be and ^B flux of the Bahcall-Ulrich SSMs (dots), and the 
90% C.L. contour of our parametrized SSM (solid curve). The magnitudes and the correlations of 
the fluxes are in excellent agreement for the two methods. 

FIG. 4. The allowed regions of the Homestake, Kamiokande, and gallium experiments, and the 
combined observations. The Earth effect is not included. The SSM uncertainties are calculated by 
(a) the Bahcall-Ulrich 1000 Monte Carlo SSMs, (b) the parametrized Bahcall-Ulrich SSM, and (c) 
the parametrized Bahcall-Pinsonneault SSM. The comparison of (a) and (b) shows that the Monte 
Carlo method and our parametrization method are essentially equivalent. The allowed regions are 
defined by ^^(sin^ 26^, Am^) < ^^.^ + A^^ with Ax^ = 6.0 at 95% C.L. 

FIG. 5. Same as Fig. ^, except that the Earth effect is included. The Kamiokande day-night 
data exclude the region shown with a dotted line (90% C.L.), and are incorporated in the combined 
fit. 
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FIG. 6. The experimental constraints on the MSW parameters assuming the Bah- 
call-Pinsonneault SSM when (a) the theoretical uncertainties of the SSM fluxes and the detector 
cross sections are ignored, (b) the theoretical uncertainties are included, but their correlations 
between the fluxes and between the experiments are ignored. Displayed in (c) is the result at 



90% C.L. calculated with the same condition in Ref |21|: the correlations and the Earth effect are 
ignored and the same experimental input parameters are used. The omission of the correlations 
overestimates the allowed large- angle regions. 

FIG. 7. The experimental constraints when the Turck-Chieze-Lopes SSM is assumed. The 
two figures are (a) without and (b) with the Earth effect. Since this SSM predicts a smaller 
flux with a larger uncertainty compared to the Bahcall-Pinsonneault SSM, the allowed regions are 
noticeably larger than those shown in Fig. ^ and ||. Especially, no constraints are obtained at 95% 
C.L. from the upper limit of the the Kamiokande result, and the Kamiokande allowed regions are 
outside the triangles with the solid lines. The region excluded by the Kamiokande II day-night 
data is shown with the dotted line (90% C.L.) in (b). 

FIG. 8. The experimental constraints (a) when theoretical uncertainties of the Bah- 
call-Pinsonneault SSM fluxes and the detector cross sections are doubled, and (b) the Homestake 
experimental uncertainty is tripled. The Earth effect and the Kamiokande II day-night data are 
included in both cases. 

FIG. 9. The Kamiokande contours of the signal to SSM ratio (including the neutral current 
reaction in flavor oscillations) for three different analytic MSW approximations of (a) Petcov, (b) 
Parke, and (c) Pizzochero. The formulae by Parke and Pizzochero fail in the large-angle, small 
Am^ region {Am? < 3 x 10^'^eV^). 
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FIG. 10. The combined allowed regions are simply taken as overlaps of the three experimental 
constraints at 95% C.L. This C.L. definition overestimates the uncertainties, allowing a parameter 
space which is marginally allowed by different experiments, but, in fact, its ^ is large and the 
combined fit is poor. As a result the obtained parameter space is significantly overestimated 
compared to Fig. |5[ especially in the large-angle region; even a third allowed region appears in 
large-angle, small Am^. Also this overlap procedure ignores uncertainty correlations between 
different experiments. 

FIG. 11. The allowed regions of the combined experiments using an improved C.L definition 
that assumes a gaussian probability density for each solution. The 90, 95, and 99% C.L. correspond 
to (a) Ax^ = 5.9, 7.3, and 10.6, ignoring the Earth effect, and (b) Ax^ = 5.5, 7.0, and 10.9, including 
the Earth effect. 

FIG. 12. The updated result of the combined MSW analysis assuming the Bah- 
call-Pinsonneault SSM. The Earth effect, Kamiokande II day-night data, theoretical uncertainties 
and their correlations are all included and the improved C.L. definition is used. 

FIG. 13. The updated result of the combined MSW analysis for the oscillations into sterile 
neutrinos. The Bahcall-Pinsonneault SSM is assumed. Bahcall-Pinsonneault SSM. The Earth ef- 
fect, Kamiokande II day-night data, theoretical uncertainties and their correlations are all included 
and the improved C.L. definition is used. 

FIG. 14. The combined allowed regions for oscillations into sterile neutrinos (a) without and 
(b) with the Earth effect. No solutions are allowed in the large-angle regions even at 99% C.L. 

FIG. 15. The experimental constraints assuming nonstandard solar models that predict neu- 
trino fluxes smaller than the SSM: (a) the opacity is reduced by 20% at temperature larger than 



5 X 10^ K m, (b) Sii is increased by 30% ^ 
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FIG. 16. The experimental constraints assuming nonstandard solar models that predict neu- 



trino fluxes larger than the SSM: (a) a high Y model |22| and (b) the maximum rate model in 



which 533 is artificially set to zero, maximizing the gallium rate to 303 SNU 

FIG. 17. The allowed regions of the combined observations when the central temperature 
(Tc) is a completely free parameter. The combined data constrain Tc to 1.02 it 0.02 (lex), which is 
consistent with the SSM (Tq = l=b0.006). Also displayed is the region excluded by the Kamiokande 
day-night data at 90% C.L.; the exclusion comes from the comparison between the different time 
bins, and is insensitive to the absolute flux or Tq- 

FIG. 18. The experimental constraints when nonstandard Su values are assumed. Sn is 
directly proportional to the flux and has the largest uncertainty among the SSM input parame- 
ters. We take 5i7 to be (a) 30% smaller than the SSM, (b) 50% larger than the SSM. In (c) the ®B 
flux is treated as a completely free parameter; the combined data constraint % = 1.43 + 0.65 — 0.42 
of the standard value (Icr) with Xmin ~ ^-^ ^^r 6 degrees of freedom. Also shown in (c) is the 
excluded region from the Kamiokande II day-night data (90% C.L.), which is independent of the 
®B flux magnitude. 

FIG. 19. (a) The combined result of the Homestake and Kamiokande experiments. From the 
allowed regions we can predict the gallium rate to be < 100 SNU at 95% C.L., which is consistent 
with the current observations of SAGE and GALLEX. (b) The combined result when the gallium 
experimental uncertainty is reduced by a factor l/\/2, which is equivalent to the data set through 
the end of 1994. The central value of the gallium rate is assumed to stay at the current value. The 
present values are used for the other experiments. 

FIG. 20. The SNO charged current spectrum expected for the nonadiabatic and large-angle 
solution. Astrophysical solutions would be similar to the large-angle spectrum. To compare the 
difference, the large-angle spectrum is normalized to the the nonadiabatic. The charged current 



cross section [142.441 and the detector resolution |39(1 are included. The error bars indicate the 



statistical uncertainties from 6,000 events (equivalent to two years of operation). 
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FIG. 21. The Super-Kamiokande electron spectrum expected for the nonadiabatic and 
large-angle solutions. To compare the difference, the large-angle spectrum is normaUzed to the 



the nonadiabatic. The detector resolution is included |40|. The error bars indicate the statistical 



uncertainties from 16,000 events (equivalent to a two year operation). 

FIG. 22. The Earth effect expected in SNO for the (a) nonadiabatic and (b) large angle 
solution. The night rate is shown with five bins according to the angle between the direction to the 
Sun and the nadir at the detector. The error bars indicate the statistical uncertainties equivalent 
to one year of operation (3,000 events). 

FIG. 23. The Earth effect expected in Super-Kamiokande for the (a) nonadiabatic solution 
and (b) large angle solution. The night rate is shown with five bins according to the angle between 
the direction to the Sun and the nadir at the detector. The error bars indicate the statistical 
uncertainties equivalent to one year of operation (8,000 events). 

FIG. 24. The Earth effect expected in Kamiokande for the (a) nonadiabatic and (b) large 
angle solution. The night time rates are shown with five bins according to the angle between 
the direction to the Sun and the nadir at the detector. The error bars indicate the statistical 
uncertainties assuming the total signals of 200 events. 

FIG. 25. The contours of the signal to the SSM ratio including the time-averaged Earth effect 
for (a) the SNO charged current mode {vg + d p + p + e), (b) Super-Kamiokande {ve scattering), 
and (c) BOREXINO {ve scattering). Superimposed are the combined current allowed regions at 
95% C.L of the Homestake, Kamiokande, and gallium experiments. 
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FIG. 26. Shown are the constraints on the Ve <-> v^i and <-> v-r oscillations |5C], the 
constraints on the oscillations into sterile neutrinos (fg <-> Vg and z^s) from big-bang nucle- 

osynthesis (BBN) [^], and three observational hints of neutrino mass: the MSW hypothesis for 
the solar neutrino deficit, the oscillation interpretation of the atmospheric neutrino deficit |46|, 
and the cold plus hot dark matter scenario for the COBE and large scale structure data |47]. Also 
displayed are predictions from various theoretical models [p!2,4^,49|. 
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FIGURES 




FIG. 12. The updated result of the combined MSW analysis assuming the Bahcall-Pinsonneault 
SSM. The Earth effect, Kamiokande II day-night data, theoretical uncertainties and their correla- 
tions are all included and the improved C.F. definition is used. 
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FIG. 13. The updated result of the combined MSW analysis for the oscillations into sterile 
neutrinos. The Bahcall-Pinsonneault SSM is assumed. Bahcall-Pinsonneault SSM. The Earth 
effect, Kamiokande II day-night data, theoretical uncertainties and their correlations are all included 
and the improved C.L. definition is used. 



50 



10^ 



10^ 



10' 



/ ' COBE 
/ j CDM+HDK/k^ 



^e-v, limit 



Atmos. 



3. 

CM 

E 




sin^2e 

FIG. 26. Shown are the constraints on the z/g ^ //^ and z/^ ^ z/^ oscillations [50], the constraints 
on the oscillations into sterile neutrinos (z/g ^ I's and z/^ ^ z/^) from big-bang nucleosynthesis 
(BBN) [35], and three observational hints of neutrino mass: the MSW hypothesis for the solar 
neutrino deficit, the oscillation interpretation of the atmospheric neutrino deficit [46], and the cold 
plus hot dark matter scenario for the COBE and large scale structure data [47]. Also displayed are 
predictions from various theoretical models [12,48,49]. 
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This figure "figl-3.png" is available in "png" format from: 
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FIGURES 




(a) pp (b) ^Be (c) 8b 

FIG. 1. The distributions of the (a) pp, (b) ''Be, and (c) ®B flux of the 1000 Monte Carlo SSMs by Bahcall-Ulrich (histograms) are flt with 
the parametrized method (solid curves) that assumes Gaussian distributions of the central temperatures and the nuclear reaction cross sections 
around their central values. 
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(a) Kamiokande (®B) (b) Chlorine (c) Gallium 

FIG. 2. The distributions of the experimental rates obtained from Bahcall-Ulrich Monte Carlo fluxes for the (a) Kamiokande, (b) chlorine, 
and (c) gallium experiments (histograms). They are compared with the rate distributions obtained from the parametrized fluxes (solid curve). 
In both cases the detector cross section uncertainties are not included. 
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FIG. 3. The distributions of the ''Be and ®B flux of the Bahcall-Ulrich SSMs (dots), and the 90% C.L. contour of our parametrized SSM 
(solid curve). The magnitudes and the correlations of the fluxes are in excellent agreement for the two methods. 
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sin^2e sin^2e sin^2e 



(a) Monte Carlo SSM (b) Bahcall-Ulrich SSM (c) Bahcall-Pinsonneault SSM 

FIG. 4. The allowed regions of the Homestake, Kamiokande, and gallium experiments, and the combined observations. The Earth effect 
is not included. The SSM uncertainties are calculated by (a) the Bahcall-Ulrich 1000 Monte Carlo SSMs, (b) the parametrized Bahcall-Ulrich 
SSM, and (c) the parametrized Bahcall-Pinsonneault SSM. The comparison of (a) and (b) shows that the Monte Carlo method and our 
parametrization method are essentially equivalent. The allowed regions are defined by x^(sin^ 2^, Am^) < Xmin "I" ^x"^ with Ax^ = 6.0 at 95% 
C.L. 
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sin^2e sin^2e sin^2e 



(a) Monte Carlo SSM (b) Bahcall-Ulrich SSM (c) Bahcall-Pinsonneault SSM 

FIG. 5. Same as Fig. 4, except that the Earth effect is included. The Kamiokande day-night data exclude the region shown with a dotted 
line (90% C.L.), and are incorporated in the combined fit. 
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(a) No theory errors (b) No correlations of theory errors (c) No correlations without the Earth effect 

FIG. 6. The experimental constraints on the MSW parameters assuming the Bahcall-Pinsonneault SSM when (a) the theoretical uncertainties 
of the SSM fluxes and the detector cross sections are ignored, (b) the theoretical uncertainties are included, but their correlations between the 
fluxes and between the experiments are ignored. Displayed in (c) is the result at 90% C.L. calculated with the same condition in Ref [21]: 
the correlations and the Earth effect are ignored and the same experimental input parameters are used. The omission of the correlations 
overestimates the allowed large-angle regions. 
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FIG. 7. The experimental constraints when the Turck-Chieze-Lopes SSM is assumed. The two figures are (a) without and (b) with the 
Earth effect. Since this SSM predicts a smaller ®B flux with a larger uncertainty compared to the Bahcall-Pinsonneault SSM, the allowed regions 
are noticeably larger than those shown in Fig. 4 and 5. Especially, no constraints are obtained at 95% C.L. from the upper limit of the the 
Kamiokande result, and the Kamiokande allowed regions are outside the triangles with the solid lines. The region excluded by the Kamiokande 
II day-night data is shown with the dotted line (90% C.L.) in (b). 
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FIG. 8. The experimental constraints (a) when theoretical uncertainties of the Bahcall-Pinsonneault SSM fluxes and the detector cross 
sections are doubled, and (b) the Homestake experimental uncertainty is tripled. The Earth effect and the Kamiokande II day-night data are 
included in both cases. 
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sin^2e sin^2e sin^2e 



(a) Petcov formula (b) Parke formula (c) Pizzochero formula 



FIG. 9. The Kamiokande contours of the signal to SSM ratio (including the neutral current reaction in flavor oscillations) for three different 
analytic MSW approximations of (a) Petcov, (b) Parke, and (c) Pizzochero. The formulae by Parke and Pizzochero fail in the large-angle, small 
Am^ region (Am^ < 3 X lO'^eV^). 
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FIG. 10. The combined allowed regions are simply taken as overlaps of the three experimental constraints at 95% C.L. This C.L. definition 
overestimates the uncertainties, allowing a parameter space which is marginally allowed by different experiments, but, in fact, its is large 
and the combined fit is poor. As a result the obtained parameter space is significantly overestimated compared to Fig. 5, especially in the 
large- angle region; even a third allowed region appears in large- angle, small Am^. Also this overlap procedure ignores uncertainty correlations 
between different experiments. 
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FIG. 11. The allowed regions of the combined experiments using an improved C.L definition that assumes a gaussian probability density 
for each solution. The 90, 95, and 99% C.L. correspond to (a) Ax^ = 5.9, 7.3, and 10.6, ignoring the Earth effect, and (b) Ax^ = 5.5, 7.0, and 
10.9, including the Earth effect. 
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sin^2e sin^2e 



(a) Sterile v without the Earth effect (b) Sterile v with the Earth effect 

FIG. 14. The combined allowed regions lor oscillations into sterile neutrinos (a) without and (b) with the Earth effect. No solutions are 
allowed in the large-angle regions even at 99% C.L. 
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(a) Low opacity model (b) High model 

FIG. 15. The experimental constraints assuming nonstandard solar models that predict neutrino fluxes smaller than the SSM: (a) the 
opacity is reduced by 20% at temperature larger than 5 X 10^ K [36], (b) 5*11 is increased by 30% [37]. 
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(a) High Y Model (b) Maximum rate model 

FIG. 16. The experimental constraints assuming nonstandard solar models that predict neutrino fluxes larger than the SSM: (a) a high Y 
model [22] and (b) the maximum rate model in which 6*33 is artiflcially set to zero, maximizing the gallium rate to 303 SNU [8]. 
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FIG. 17. The allowed regions of the combined observations when the central temperature (Tc) is a completely free parameter. The combined 
data constrain Tc to 1.02±0.02 (It), which is consistent with the SSM (Tc = 1±0.006). Also displayed is the region excluded by the Kamiokande 
day-night data at 90% C.L.; the exclusion comes from the comparison between the different time bins, and is insensitive to the absolute ®B flux 
or Tc. 
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FIG. 18. The experimental constraints when nonstandard 5*17 values are assumed. 5*17 is directly proportional to the ®B flux and has the 
largest uncertainty among the SSM input parameters. We take 5*17 to be (a) 30% smaller than the SSM, (b) 50% larger than the SSM. In (c) 
the ®B flux is treated as a completely free parameter; the combined data constraint = 1.43 + 0.65 — 0.42 of the standard value (It) with 
X^i^ = 2.0 for 6 degrees of freedom. Also shown in (c) is the excluded region from the Kamiokande II day-night data (90% C.L.), which is 
independent of the ®B flux magnitude. 
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(a) Kam & CI combined (b) Ga projection (1994) 

FIG. 19. (a) The combined result of the Homestake and Kamiokande experiments. From the allowed regions we can predict the gallium 
rate to be < 100 SNU at 95% C.L., which is consistent with the current observations of SAGE and GALLEX. (b) The combined result when 
the gallium experimental uncertainty is reduced by a factor 1/^2, which is equivalent to the data set through the end of 1994. The central 
value of the gallium rate is assumed to stay at the current value. The present values are used for the other experiments. 
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FIG. 20. The SNO charged current spectrum expected for the nonadiabatic and large-angle solution. Astrophysical solutions would be 
similar to the large-angle spectrum. To compare the difference, the large-angle spectrum is normalized to the the nonadiabatic. The charged 
current cross section [43,44] and the detector resolution [39] are included. The error bars indicate the statistical uncertainties from 6,000 events 
(equivalent to two years of operation). 



57 



Super-Kamiokande Electron Spectrum 
1 6000 Events (2 year) 




5 10 15 

E3 (MeV) 



FIG. 21. The Super-Kamiokande electron spectrum expected for the nonadiabatic and large-angle solutions. To compare the difference, the 
large- angle spectrum is normalized to the the nonadiabatic. The detector resolution is included [40]. The error bars indicate the statistical 
uncertainties from 16,000 events (equivalent to a two year operation). 
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FIG. 22. The Earth effect expected in SNO for the (a) nonadiabatic and (b) large angle solution. The night rate is shown with five bins 
according to the angle between the direction to the Sun and the nadir at the detector. The error bars indicate the statistical uncertainties 
equivalent to one year of operation (3,000 events). 
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FIG. 23. The Earth effect expected in Super-Kamiokande for the (a) nonadiabatic solution and (b) large angle solution. The night rate 
is shown with five bins according to the angle between the direction to the Sun and the nadir at the detector. The error bars indicate the 
statistical uncertainties equivalent to one year of operation (8,000 events). 
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FIG. 24. The Earth effect expected in Kamiokande for the (a) nonadiabatic and (b) large angle solution. The night time rates are shown 
with five bins according to the angle between the direction to the Sun and the nadir at the detector. The error bars indicate the statistical 
uncertainties assuming the total signals of 200 events. 
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(a) SNO (charged current) (b) Super-Kamiokande (c) BOREXINO 

FIG. 25. The contours of the signal to the SSM ratio including the time- averaged Earth effect for (a) the SNO charged current mode 
(z/e + d ^ p + p + e), (b) Super-Kamiokande (z/e scattering), and (c) BOREXINO (//e scattering). Superimposed are the combined current 
allowed regions at 95% C.L of the Homestake, Kamiokande, and gallium experiments. 
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